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Abstract 



We are interested in the long-time asymptotic behavior of growth-fragmentation equations with 
' a nonhnear growth term. We present examples for which we can prove either the convergence to a 

' steady state or conversely the existence of periodic solutions. Thanks the General Relative Entropy 

method applied to well chosen self-similar solutions, we show that the equation can "asymptotically" 
be reduced to a system of ODEs. Then stability results are proved by using a Lyapunov functional, 
and existence of periodic solutions are proved thanks to the Poincare-Bendixon theorem or by Hopf 
bifurcation. 
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1 Introduction 

We are interested in growth models which take the form of a mass preserving fragmentation equation 
complemented with a transport term. Such models are used to describe the evolution of a population in 
which each individual grows and splits or divides. The individuals can be for instance cells [4} [5| 134^ |45] 
or polymers [TJ [T71 [31] and are structured by a variable x > which may be size [211 [22], label 
[2], protein content [El [42], proliferating parasites content [3] ; etc. More precisely, we denote by 
u = u{t, x) > the density of individuals of structured variable x at time t, and we consider that the 
time dynamic of the population is given by the following equation 

■§iu{t, x) + -§^{T{t, x)u{t, x)) + fi{t, x)u{t, x) = {Fu){t, x), t > 0, x > 0, 

n(t,2; = 0) = 0, i>0, (1) 
u{t = 0, x) = uq{x) > 0, rc > 0, 

where is a mass conservative fragmentation operator 

j'OO 

{Tu){t,x)=l b{t,y,x)u{t,y)dy - /3{t,x)u{t,x). (2) 
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The mass conservation for the fragmentation operator requires the relation 

P{t,x) = I ^b{t,x,y)dy. (3) 
Jo X 

The coefficient y) > represents the rate of splitting for a particule of size y at time t and 
y,x) > represents the formation rate of a particule of size x < y after the fragmentation. The 
velocity r(t, x) > in the transport term represents the growth rate of each individual, and fi{t, x) > 
is a degradation or death term. 

We consider that the time dependency of r and is of the form 

T{t,x) = V{t)T{x) and fi{t,x) = R{t)fi{x) (4) 

and moreover that the size dependency is a powerlaw 

r(x) = Tx'^ and fJ-ix) = fi. (5) 

The choice of the coefficients V{t) and R{t) depends on the cases we want to analyse. We give below 
four examples in which they are nonlinear terms or periodic controls. The fragmentation coefficients 
are assumed to be time-independent and of self-similar form 

Kt,x,y) = ^.(l) (6) 
X \xJ 

with K a nonnegative measure on [0, 1]. For h under this form, the quantity 

no := / k{z) dz 
Jo 

represents the mean number of fragments produced by the fragmentation of an individual. We remark 
that under Assumption ([6]), relation ([3]) becomes z k{z) dz = 1 and then it enforces uq > 1. If k 
is symmetric {k{z) = k(1 — z)), we have necessarily uq = 2. We also assume that /3 is a powerlaw 
coefficient 

P{x) = I3x^ (7) 

and we denote by the fragmentation operator associated to coefficients satisfying Assumptions ([6|) 
and dZD. 

Now we state our main results concerning different choices for V{t) and R{t). First we investigate the 
nonlinear growth- fragmentation equations corresponding to the case when V or /and R are functions 
of the solution u{t, x) itself. We also consider a model of polymerization in which the transport term 
depends on u and on a solution to an ODE coupled to the growth-fragmentation equation. The long- 
time asymptotic behavior of these equations is investigated under the assumption that t(x) is linear 
(i.e. u = 1) and (3 increasing (i.e. 7 > 0). We finish with a study of the long-time asymptotics in the 
case when V and R are known periodic controls. 
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Example 1. Nonlinear drift-term. We consider that the death rate is time independent {R = 1) 
and that the transport term depends on the solution itself through the nonlinearity 

d / f \ d 

— x) = —f i / x^u{t, x) dx j — (x u{t, x)) — fi u(t, x) + Fr^u{t, x) (8) 

where / : — t- is a continuous function which represents the influence of the wheighted total 
population J x^u(t, x) dx [p > 0) on the growth process. Such weak nonlinearities are common in 
structured populations (see [361 EH EO] for instance). The stability of the steady states for related 
models has already been investigated (see [iHl [261 [271 [281 [50]) but never for the growth- fragmentation 
equation with the nonlinearities considered here. We prove in Section [3] convergence and nonlinear 
stability results for Equation ([8]) in the functional space H := L'^{{x + x^)dx) for r large enough, 
or more precisely in its positive cone denoted by 7^"*". These results are stated in the two following 
theorems. They require that / is continuous and satisfies 

M := {/; /(/) = /u} is a finite set and limsup/(/) < fi. (9) 

Theorem 1.1 (Convergence). Assume that f satisfies Assumption ([9|), that 7 G (0,2] and that the 
fragmentation kernel k satisfies Assumption (|32p . Then the number of positive steady states of Equa- 
tion ([8|) is equal to jjA/" and any solution with an initial distribution uq € Ti^ either converges to one 
of these steady states or vanishes. 

Theorem 1.2 (Local stability). Assume that f G C^(M4.) satisfies Assumption ([9|), that 7 S (0,2] 
and that the fragmentation kernel k satisfies Assumption (j32p . Then the trivial steady state is locally 
exponentially stable if f{0) < fJ- and p > 1, and unstable if f{0) > /i. Any nontrivial steady state Uoo is 
locally asymptotically stable if f (/ xPuoo{x) dx) < 0, locally exponentially stable if additionally k = 2, 
and unstable if f (J x^Uoo{x) dx) > 0. 

As an immediate consequence of these two theorems, we have the following corollary. 

Corollary 1.3 (Global stability). For 7 G (0,2] and under Assumption (f32]l . if f € C^(M+) sat- 
isfies Q with N a singleton, then there is a unique nontrivial steady state Uoo, and if additionally 
f (/ x^Uoo{x) dx) < 0, then it is globally asymptotically stable in ^"*"\{0}. 

The method combines several arguments. First it uses the General Relative Entropy principle 
introduced by |48[ l49l [51] for the linear case. Secondly it reduces the system to a set of ODEs which 
has the same asymptotic behavior than Equation ([8]). Then we build a Lyapunov functional for this 
reduced system. Therefore our result extends to the case of the growth-fragmentation equation several 
satability results proved for the nonlinear renewal equation in \A7\ I53j. 

Example 2. Nonlinear drift and death terms. We can also treat several nonlinearities as in 

—u{t,x) = —f^J x^u{t,x) dx^ —{^xu(t,x)) — g (^j x'^u{t,x) dx^ u{t,x) -\- F^u{t,x). (10) 

In this case, we show that oscillating behaviors can appear. The existence of non-trivial periodic 
solution for structured population models is a very interesting and difficult problem. It has been 
mainly investigated for age structured models with nonlinear renewal and/or death terms, but there 
are very few results [H [6l [IH [381 HHl HH IHl 153 [58]. Taking advantage of the Poincare-Bendixon 
theorem, we prove the existence of oscillating solutions for Equation (jlOp and the "stability" of this 
behavior in a sense specified in the following theorem. 
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Theorem 1.4. Forj G (0,2] and under Assumption ([32]) . there exist functions f and g and parameters 
p and q for which we can find an open set V C Ti^ with the property that any solution to Equation (jlOh 
with an initial distribution uq £ V converges to a periodic solution. 

The proof of this resuh is given in Section [H 



Example 3. The prion equation. In Section [5l we are interested in a general so-cahed prion 
equation 

dV{t) 



dt 



j xPy^ j xu{t, x) dx - 5V{t) + A , 



-V{t)f (^j x^u^ (xn(t, x)) — fiu{t, x) + F^u{t, x). 



(11) 



In this equation, the growth term depend on the quantity V{t) of another population (monomers for 
the prion proliferation model). We prove for this system the existence of periodic solutions under some 
conditions on /. In age structured models, such solutions are usually built using bifurcation theory, 
particularly by Hopf bifurcation (see |41j for a general theorem). Here we consider the power p as, a. 
bifurcation parameter and we prove existence of periodic solution by Hopf bifurcation. 

Theorem 1.5. There exist a function f and parameters for which Equation (llip admits periodic 
solutions. 



Example 4. Perron vs. Floquet. Our method in Section [2.21 can also be applied to the situation 
when V{t) and R{t) are periodic controls 

d d 

Truit, x) + V(t) — (TX u(t, x)) + R(t)nu(t, x) = F^uit, x). (12) 
ot ox 

In this case. Theorem 12.11 allows to built a particular solution called the Floquet eigenvector, starting 
from the Perron eigenvector which correspond to constant parameters. Moreover, we can compare the 
associated Floquet eigenvalue to the Perron eigenvalue. The results about this problem are stated in 
Section [H 

Before treating the different examples, we explain in Section [2] the general method used to tackle 
these problems. It is based on the main result in Theorem 12 . 1 1 and the use of the eigenelements of the 
growth-fragmentation operator together with General Relative Entropy techniques. In Sections [3l HI [5] 
and[6l we give the proofs of the results in Examples 1, 2, 3 and 4 respectively. 



2 Technical Tools and General Method 
2.1 Main Theorem 

The proofs of the main theorems of this paper are based on the following result which requires to 
consider that t{x) is linear [i.e. z/ = 1). In this case, there exists a relation between a solution to 
Equation ([TJ with time-dependent parameters (Vi(t), Ri{t)) and a solution to the same equation with 
parameters {V2{t) , R2{t)) ■ More precisely, we can obtain one from the other by an appropriate dilation. 
The following theorem generalizes the change of variable used in [25] to build self-similar solution to 
the pure fragmentation equation. 



4 



Theorem 2.1. Consider that Assumptions dH)-© are satisfied with v = 1 and 7 > 0. For ui{t,x) a 
solution to Equation ([T|) with parameters {Vi,Ri), the function U2{t,x) defined by 

U2{t,x) = W-''{t)ui (^h{t),W-''{t)x^ ^^.J^{W{s)R^{syR2is))ds ^^3) 

with k = ^, is a solution to Equation ^ with {V2,R2) ifW>0 and h satisfy 

W = TK^V2-ViW), (14) 
h = W. 

Conversely, if h : M_|_ — ^ M4. is one to one and if U2 is a solution with (V2, R2), then ui defined by (|13|) 
is a solution with (Vi,iii). 

The proof of this result is nothing but easy calculation that we leave to the reader. 



Remark 2.2. To check that h is one to one, we can take advantage that ODE (jl4p is a Bernoulli 
equation which can be integrated in 

^ l + WolfiV,{s)elfoy^(^')'^^'ds' ^^^^ 

To tackle the different examples, we use Theorem 12.11 together with two techniques appropriate for 
this type of equations. First we recall the existence of particular solutions to the growth-fragmentation 
equation in the case of time-independent coefficients. They correspond to eigenvectors to the growth- 
fragmentation operator and we can give their self-similar dependency on parameters in the case of 
powerlaw coefficients (see [30^ [9]). This dependency is the starting point which leads to Theorem 12.11 
and it allows to built an invariant manifold for Equation ([1]) in the case = 1. It also provides 
interesting properties on the moments of the solutions when u ^ 1. Then we recall results about the 
General Relative Entropy (GRE) introduced by [481 B9l [5T] for the growth-fragmentation model. This 
method ensures that the particular solutions built from eigenvectors are attractive for suitable norms. 

2.2 Eigenvectors and Self-similarity: Existence of an Invariant Manifold 

When the coefficients of Equation ([T]) do not depend on time, one can build solutions {t, x) 1— )■ U{x)e^^ 
by solving the Perron eigenvalue problem 

AU{x) = --§^{t{x)U{x)) - fiix)Kix) + iJ^U)ix), x ^ 0, 

(16) 

tU{x = 0) = 0, U{x) > 0, U{x)dx = 1. 

The existence of such elements A and lA has been first studied by [52] and is proved for general 
coefficients in [23]. The dependency of these elements on parameters is of interest to investigate the 
existence of steady states for nonlinear problems (see |12[ [9]). In the case of powerlaw coefficients, 
we can work out this dependency on frozen transport parameter V and death parameter R (see [30]). 
Under Assumptions (I5])-(I7|), the necessary condition which appears in \12>\ I46j to ensure the existence 
of eigenelements is 7 + 1 — > 0. Then we define a dilation parameter 

k := > (17) 

7 + 1-1/ 
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and we have explicit self-similar dependencies 



A{V, R) = V''^A{1, 0) - Rfi and U{V; x) = V'^Uil; V'^x). (18) 

The eigenvector U does not depend on i?, that is why we do not label it. Thereafter A(1,0) and 
•) are denoted by A and U for the sake of clarity. The result of Theorem 12. II is based on the idea 
to use these dependencies to tackle time-dependent parameters. An intermediate result between the 
formula (jl3p and the dependencies (jlSp is given by the following corollary. 



Corollary 2.3. Under the assumptions of Theorem \2. 11 if W is a solution to 

W=^(V-W) (19) 

k 

then 

u{t, x) = U{W{ty, x)e^o A{W{s),R{s))ds (20) 

is a solution to Equation ([T|). 

Proof. For = 1, we can compute A = A(1,0) by integrating Equation (116p with ^ = against xdx. 
We obtain thanks to the mass preservation of 



fOO /'OO 

A / xU{x)dx = T xU{x)dx 
Jo Jo 



and so A(l, 0) = r. Thus using the dependency (fTSj) we find that A(W, 0) = tW and Equation (fT9]) is 
nothing but a rewriting of Equation ()14p . We use this formulation (I19p here to highlight the link with 
eigenelements, and because it allows after to obtain results in the cases when ^ 1. Now we apply 
Theorem 12.11 for Vi = 1, Ri = and V2 = V, R2 = R and we obtain that 

U2{t,x) = W-''U{W-''x)e^^^^^-^oRi^)d' =u{W{ty,x)e-foMWis),Ris))ds 
is a solution to Equation ([T|). □ 



This corollary provides a very intuitive explicit solution in the spirit of dependencies (jlSp . At each 
time t, the solution is an eigenvector associated to a parameter W{t) with an instantaneous fitness 
A{W{t),R(t)) associated to the same parameter W{t). The function t 1— t- W{t) thus defined follows 
V{t) with a delay explicitely given by ODE (fT9|) . 

A very useful consequence for the different applications is the existence of an invariant manifold for 
the growth- fragmentation equation with time-dependent parameters of the form ([1]). Let define the 
eigenmanifold 

£ ■.= {QU{W;-), {W,Q)e{M.lf} (21) 

which is the union of all the positive eigenlines associated to a transport parameter W. Then Corol- 
lary [23] ensures that, under the assumptions of Theorem 12.11 any solution to Equation ([TJ with an 
initial distribution uq £ £ remains in £ for all time. Moreover the dynamics of such a solution reduces 
to the ODE (jl9p and this is the key point we use to tackle nonlinear problems. 



For 7^ 1, the technique fails and we cannot obtain explicit solution with the method of Theorem 12. 11 
Nevertheless, we can still define W as the solution to ODE (jl9p and give properties of the functions 
defined as dilations of the eigenvector by ()20p . We obtain that the moments of these functions 
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satisfy equations that are similar to the ones verified by the moments of the solution to the growth- 
fragmentation equation. In the special case u = 0, and 7 = 1, we even obtain the same equation. 
More precisely, if we denote, for a > 0, 

Ma[W]it):= x''UiWit);x)e^oMWis)Ms))ds^^^ ^22) 
Jo 

and 



MJu] 



i]{t) := / x'^u{t,x)dx, (23) 
then we have the following result. 

Lemma 2.4. On the one hand, if W is a solution to Equation (I19p . then the moments Ma satisfy 

Ma = akOaVMa+u-l + (1 " a)hbaMa-^-i - fiRMa (24) 

with 

Ma[U] , , Ma[U] 

O-a '■= T7 FTTT ^""^ "a ■~ 



On the other hand, if u is a solution to the frag entation- drift equation, then the moments Ma satisfy 

M„ = arVMa+u-i + (ca - l)/3Ma+-y - iiRMa (25) 
with ^ 

Ca '■= / z°'k{z) dz. 



JO 

Proof. Thanks to a change of variable we can compute for all a > 

Ma[W] = Ma[U]W''"e-fo^^^('^'R('^^ds 

so we have 

W 

Ma = ka—Ma + MW,R)Ma 
W 

= aAW^^-\V - W)Ma + AW'^^Ma - f^RMa 
= aAVW^'^'^-^^Ma + (1 - a)AW''''Ma - f^RMa 

= aAVOaMa+u-l + (1 - a)AbaMa+'y - llRMa- 

Integrating Equation ([1]) against x^dx we obtain by integration by part and thanks to the Fubini 
theorem 

d 



dt 



x°'u{t,x)dx = tV / x'^dx{x'^u{t,x)) dx — fiR / x°'u{t,x)dx 



CO fOO 



-f3 I x°+Tu(t, x)dx + p I x" / y^-^K ( - ) dydx 



arV J x'^+''-'u{t,x)dx - fiR J x°'u{t,x)dx 

-/? / x"+^u{t,x)dx + /3 ry^+^ r^'^(-) —dy 



'0 y \vj y 

arVMa+u-i + {Ca - l)PMa+j - flRMa- 



□ 
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In the particular case = 0, 7 = 1 and k symmetric, we can compute the constants Oc, ha and Cq, 
for a = 1 or a = 2. A consequence is the useful result given in the following corollary. 

Corollary 2.5. In the case when = 0, 7 = 1 and n symmetric, the zero and first moments 
and (Mo,Mi) are both solution to 



U 
P 



tV -i^iR 



U 
P 



(26) 



This Corollary allows in Section [6] to compare the Perron and Floquet eigenvalues not only for u = 1 
but also for = 0, 7 = 1. We do not have in this case a particular solution to the growth- fragmentation 
equation as in Corollary 12.31 but a particular solution of the reduced ODE system (126p . 

Proof. For k symmetric, we have already seen that cq = hq = 2. Together with Assumption ([3]) which 
gives ci = 1, we conclude that {Mq,Mi) is solution to Equation (j26l) . 
Integrating Equation (I16p against dx and xdx we obtain 



It allows to compute 



A = /3 J xU{x)dx and A j xU{x)dx = T. 
j xU{x) dx = and A = \/r^. 



Thus ai = . /5 and 60 = yf, and (A^Oi-^i) satisfies Equation ()26p thanks to Lemma [27 



□ 



2.3 General Relative Entropy: Attractivity of the Invariant Manifold 

The existence of the invariant manifold £ is useful to obtain particular solutions to nonlinear growth- 
fragmentation equations. But what happens when the initial distribution uq does not belong to this 
manifold ? The GRE method ensures that £ is attractive in a sense to be defined. 

The GRE method requires to consider the adjoint growth- fragmentation equation 



d d 

— \l}{t,x) = T{t,x)—i){t,x) - ll{t,x)\l){t,x) + {J^*i)){t,x) 



where J-* is the adjoint fragmentation operator 



(27) 



(^»(t, x) := / 6(t, X, y)i^{t, y) dy - /3(t, x)ij{t, x). 
Jo 

If u and V are two solutions to Equation ([1]) and ^ is a solution to Equation ()27p . then we have for 
any function if : M — > M 



00 roo 



Jy 



b{t,y,x)ip{t,x)v{t,y) 



H 



u{t, x) 
v{t, x) 



ir/^ii^y)\ I ( u{t,x)\ ( u{t,y) u{t,x) 



v{t,x)J \v{t,y) v{t,x) 



dxdy. 
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When H is convex, the right hand side is nonpositive and we obtain a nonincreasing quantity called 
GRE. 



In the case of time-independent coefficients, we can chose for v a solution of the form U{x)e^* . Then, 
to apply the GRE method, we need a solution to the adjoint equation and such solutions are given by 
solving the adjoint Perron eigenvalue problem 

K<j){x) = T{x)^{<i){x)) - ii{x)c^{x) + {F*ct>){x), X ^ 0, 

(28) 

<P{x) > 0, (t){x)U{x)dx = 1. 

Such a problem is usually solved together with the direct problem (|16|) and the first eigenvalue A 
is the same for the two ones (see |231 HGj 151]). Then the GRE ensures that any solution u to the 
growth-fragmentation equation behaves asymptotically like lA{x)e^^. More precisely it is proved in 
|491 [51] under general assumptions that 

^hm \\Q^Mtr)e-^' -U\\Lv{W'P4>dx) =0, Vp > 1. (29) 

where Qo = J uo{y)(l){y) dy with no(x) = u(t = 0,x). 

Under the assumptions of Theorem 12.11 and for p = 1, this convergence result can be interpreted as 
the attractivity of the invariant manifold £ in L^{(j)dx) with the distance 

d{u,£) := inf \\q~^u - ^||Li(^(a;)dx) 

where Q ■= J u{y)4>{y) dy. Consider, for V{t) > 0, a solution W to 

W 

W = ^(t-W) 
k 

with W{Q) = 1. First we have W > -^W"^, so W > and h > kln{l + ^). Thus /i : M+ M+ is 
one to one and we can define from a solution u{t, x) the function 

v{h{t),x) := H^'=(t)n(t,W^'=(t)x)e-i'o(W'W-/^«W)rf'* (30) 

which satisfies v{0,x) = u{0,x) = uo{x). Then, using Theorem 12.11 v{t,x) is a solution to 

d d 

—v{t,x) = - — {xv{t,x)) - v{t,x) + {J^^v){t,x) 

so we have 

foo poo 

v{t,x)(j){x) dx = / v{t = 0,x)(f){x) dx = Q(j 



and 

lim ||^o^t'(t,-) -^llLi(</,(x)dx) = 0. 
For = 1, is linear (see examples in [23]), so we can compute from ([30]) 



g{t)= [ u{t,ymy)dy = goW''{t)eIoims)-,Ris))ds 
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and 

d{u{t,-),£) < \\Q-'^{t)u{t,-) - U\\ = Wg-^^vit,-) - U\\ 0. (31) 

The exponential decay in (I29p is proved in |52l [39] for p = 1 and for a constant fragmentation rate 
/3(x) = /3. It is also proved in [8] for powerlaw parameters in the norm corresponding to p = 2 and 
this is the case we are interested in. Assume that the coefficients satisfy ([5]), ([6]) and ([7]) and assmne 
also that the fragmentation kernel is upper and lower bounded 

3k,k>0, VzG[0, 1], k<k{z)<k. (32) 

Then a spectral gap result is proved in [8] for u = 1 in L^(U~^(j)dx) and then the result is extended 
to bigger spaces thanks to a general method for spectral gap. 

Theorem [8]. Under Assumption ([5]) with u = 1, Assumptions ([6]) and (j32p . and Assumption ([7|) 
with 7 G (0, 2], there exists f > 3 and for any a G (0,7) and any r > r there exists Ca^r such that for 
any uq £ Ti := L"^ (6) , 6{x) = x + x'', there holds 

Vt>0, \\Qo^uit,-)e-^' -U\\n<Ca,r\\6o^uo-U\\ne-'''. (33) 

This result is very useful for Applications 1 and 2 because L^{9) C L^{xP) for r > 2p + 1. Moreover 
the exponential decay allows to prove exponential stability results for Equation ([8]) when k is constant 
(see Section [3|). 



3 Nonlinear Drift Term: Convergence and Stability 

Consider the nonlinear growth-fragmentation equation ([8]) where the transport term depends on the 
p-moment of the solution itself. This dependency may represent the influence of the total population 
of individuals on the growth process of each individual. We study the long-time asymptotic behavior 
of the solutions in the positive cone with the weight 6{x) = x + x^ for 

r > max(f, 2p+ 1). (34) 

We prove that there is always convergence to a steady state, provided that the function / is less than fi 
at the infinity. This result is precised in Theorem 13.11 which immediately leads to Theorem 11.11 and the 
stability of the steady states is given in Theorem ll.2l We use the notation Mp for Mp[l/(] = J x^ZY(x) dx. 

Theorem 3.1. Assume that f is a continuous function on [0,+oo) which satisfies Assumption ([9]), 
that 7 G (0,2] and that the fragmentation kernel k satisfies Assumption (I32p . Then the nontrivial 
steady states of Equation ([8]) write 



^/i-^P^/(/i;-) with loo^Af, 
and for any solution u, there exists loo G A/" U {0} such that 

= 0. (35) 



lim 



^(t,-)-^/^"'^^(/x; 



H 
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Proof of Theorem \3.1[ First step: uq G £. 

Consider an initial distribution uq £ £ defined in ()2ip . then there exist Wq > and Qo >0 such that 



uo{x) = QoU{Wofi;x). 
Let u{t, x) be the solution to Equation ^ and define W as the solution to 

W = ^ ( f ( [ xPu{t,x)dx] - fiW 



k 

W{0) = Wq. 
Then Corollary ()2.3p ensures that 

yt,x> 0, u{t,x) = Qoli{W{t)fi;x)e^'^o(W(')-^)'i' 

and so we have 



(36) 



r xPu{t,x)dx = QoW^V^ (^j\pU{x)dx^ gM/oW^)-i)d._ 



Now defining 

we obtain a system of ODEs equivalent to Equation (l8|) in 

W 



W = -(^fp(w'^PQ)-^,w), 
Q = f,Q{W-l), 



(37) 



with the notation 

xPU{x) dx ] , 



and proving convergence of u is equivalent to prove convergence of (Wi Q). To study System ([57]) . we 
change the unknown to Z := W^^Q. Then {W, Z) is solution to 



W = —{fp{Z)-^LW), 



w 

k " (38) 

Z = pZ{fp{Z)- tiW) + f,Z{W-l), 



and the positive steady states satisfy Woo = ^/p(^oo) = 1. If p = 1, ^ is solution to Z = Z(/i(Z) — /i) 



SO for Zq > 0, Z converges to Z^o = k^,^ with /qo € A/" U {0}. When Z^o > 0, then W ^ 1 because 



fi{Z) ^ /i. 

In the case when p ^ 1, we write System (j38p as 



= -lW{f,-f,{Z))-!^,WiW-l) ^^^^ 
Z = -pZ{fi-fp{Z))-{p-l)fxZ{W-l), 
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and we exhibit a Lyapunov functional for this equation. We define 

G{W) ■=W -1- IniW) 

and 



F{Z):= [ in-Uiz))-, 



then we look for a parameter a such that k^G{W) + a^F^Z) is a Lyapunov functional for Equa- 
tion We have 



^ {k^lG{W) + a^F{Z)) = -f,\W - 1)2 - aV/i - fp{Z)f - ^(1 + a\p - 1))(1^ - l)(/i - 
so, if we find a such that 1 + a'^{p — 1) = 2ay/p, we will have 

^ (/./iCW + a^FiZ)) = -{ti{W - 1) + aVP(/^ - fp{Z))f < 0. 

For p ^ 1, there are two roots to the binomial 

{p - l)a'^ - 2y/^a + 1 = 

which are 

V^±l 1 

a — — 



p-1 ^±1 

For these two values, the dissipation is nonpositive but we want to have negativity outside of the 
steady states in order to conclude to convergence. We use a combination of this two values and the 
following lemma. 

Lemma 3.2. For any a ^ 1, there exists lo > such that 

ya,b, {a + bf + {a + abf yuia"^ + b'^). (40) 
Proof. We prove the inequality 

(a + bf + (a + abf > f{a - l){a'^ + b^) (41) 

where 

, . _ 4 - 2V4 + + 

2 + f^/4T^+f 

This function is positive except at rr = where /(O) = 0. Its maximum is /(— 2) = 2. 
To prove (jH]), we define /3 := f and we prove that 



/(I — a) = min 



(/3 + 1)2 + (/3 + a)^ 



/3 V /32 + I 

Define 

_ (/j + l)2 + (/3 + a)2 
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compute 

9a[P) (l + /32)2 

and we find that the minimum of g is reached for 

/3a = (a - 1 + V4+(a-l)2) . 
As a consequence, for ah a G M and 6 7^ 0, we have 

(a + bf + (a + a6)2 > ga{/3a){a^ + 6^) 
and QaiPa) = /(I — o)- This is stih true for 6 = by passing to the hmit. □ 

Denoting a+ = and a_ = we define L{W,Z) := 2kfj,G{W) + (a2_ + a'i)F{Z). Then 

Lemma 13.21 ensures the existence of a; > such that 



^^L{Wit),Z{t)) = -{fiiW - 1) + a+^if, - aZ))f - {^i{W - 1) + a„VP(/^ - 

< _^ (^2(^ _ -^)2 ^ ^^^(^ _ /^(^))2) z), (42) 

and -D(VF, is positive outside of the steady states. Moreover Assumption Q ensures that L and 
D are coercive in the sense that L{W,Z) — ;> +00 and D{W,Z) — )• +00 when Z)|| — ;> +00. So L 
is a Lyapunov functional for System (|39p and we can conclude to the convergence of the solution to 
a steady state. If /(O) > ^, L(VF, Z) +00 \iW ov Z tends to 0, so for any {Wq > 0, Zq > 0) the 
solution {W,Z) converges to a critical point of L, namely (1, jjiff) with /oo G AA. If /(O) < ;U, then 

for any 1^ > we have that L{W, Z) ^ — —00. So either [W, Z) converges to (1, jj^ff) with 
/oo £ AA, or Z — 0. To conclude to the convergence mTi, we write 



00 



\\u{t,-)- ZooU{^r,-)\\= / {Q{t)U{W{t)fi;x)- ZocU{fi;x)y{x + x'')dx 

Jo 

and we use dominated convergence. We know by Theorem 1 in [23] that under Assumption ()32p and 
for = 1, x"hl{x) is bounded in IR+ for all a > 0, so it ensures that the integrand is dominated 
by an integrable function. Then the ponctual convergence is given by the convergence of {W,Z), so 
convergence (f35]l occurs. 

Second step: general initial distribution uq. 

Departing from u a solution to this Equation ([8]) not necessarily in if^, we define as in Section 

v{h{t),x) := l^(t)^''u(t,t^(t)'=x)e^(*"''W) 

with 



and 

h = W. 
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We recall that in this case h : 1R+ — )• is one to one. We choose the initial values VF(0) = 1 and 
h{0) = to have v{0,x) = u{0,x) = uo{x). Then v{t,x) is a solution to 



dt dx 
and, thanks to the GRE, we conclude that 



d d 
v{t, x) = —fJ--^ [xv(t, x)) — fiv{t, x) + J^yV{t, x) 



(43) 



v{t,x) qqU{^;x) 



i— S-oo 



where 



and so 



£•0 



(t){^] x)vq{x) dx = I x)uq{x) dx 



j xPu{t,x)dx go/^ (^j xPU{x)dx^ T^(t)'=Pe^(''W-*)dx. 
As a consequence it holds 



W{t) 



w 



fp 



with Q{t) = qqc^^^^^^ *^ which satisfies the equation 

Q = fiQ{W-l). 



The interpretation of this is that System (j36p represents asymptotically the dynamics of the solutions 
to Equation ([8]). More rigorously, define 



Then we have 



e{t) :- 



W 



f xPu{t, x) dx 
fi>'PMpQ{t)W''P{t) 



1. 



(44) 



Wit) = — (/p (w'^PQil + e{t))J - fiW 



and, using the Cauchy-Schwarz inequality and the exponential decay theorem of [8], 

^-A:pg/i(t-h(t)) 



1 



< 



< 



The function is integrable under Assumption (I34p . Finally the dynamics of the solution u to 

Equation ([8]) is given by 



1 

fi'^PAdp 

Cp 
n'^PMp 



QQ^it, x) - U{fi; x)T^'=fe'^('^W-*))xP dx 
{QQ^v{h{t),x) —U{ij;x))x^ dx 

QQ^v{h{t),x) — U{ij; x)p((/>(x) + x^) dx 



,2p 



cf){x) + x' 



-dx 



(45) 



W = -lW{fi-fp{{l+e)Z))-!^W{W-l) 
Z = -pZ{fi-fp{{l + e)Z))-{p-l)fiZ{W-l) 



(46) 
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where Z = W'^pQ, and e{t) > 0. Now we look what becomes the Lyapunov functional of the first 

step for this system and we obtain 

^L{W, Z) < -DiW, Z) + {^i{W - 1) + {al + - f^{Z))) (/^((l + e)Z) - fp{Z)) . (47) 
at V ^ , 

■.=E{W,Z,e) 

Thanks to Assumption we know that fp is bounded, so W which is solution to 

W = ^{fp{{l + e)Z)-f^W) 

is also bounded. Thus E{W, Z, e) is bounded and, because L and D are coercive, the trajectory {W, Z) 
is bounded. Moreover E{W{t), Z{t), e{t)) > because e{t) — )• 0, so {W,Z) converges to a steady 

state (l,Zoo). Finally we write 

\\u{t, •) - ZMf^; Oil < Mt, •) - Q{t)U{W{t)fi; oil + \\Q{t)U{W{t)fi; •) - Z^Uifi; ■)\\ 

and we know from the first step that \\Q(t)U{W{t)ii; •) — Zooh({n; •)|| — )• 0. We treat the other term 
thanks to the spectral gap theorem of [8] 

\\u{t,-)-Qm{W{t)fx;-)\\n = qI^I \Qo'v{Ht),x) -U{fi;x)\\x + W^'-'^''x')dx^ 

< C\\g^MHt),-)-ll{f^;-)\\ 

< C\\0,\o-U{f^r)\\e-''''^'l (48) 
So ||n(t, •) — Q{t)U{W{t)fi; •)|| — )• because h{t) — )• +oo, and the proof is complete. □ 



Proof of Theorem \1.2l Stability of the trivial steady state. 

We start with the stability of the zero steady state when /(O) < fi and p > 1. For this we integrate 
Equation ([8]) against x'^ and find 



d 



xPu{t, x) dx \ < i f [ / xPu{t, x) dx] — fi] [ / xPu{t, x) dx 



thanks to the mass conservation assumption ([3|). Thus / xPu{t,x) dx is a decreasing function if 
/ (Mp[iio]) < /U. Then we integrate against x + x^ for any r > 1 and we obtain 



u{t, x) {x + x^')dx I < ( / ( / x^M(t = 0, x) I — I ( / u{t, x) {x + x^)dx 



d_ 
It 

which ensures the exponential convergence 

\\u{t,-)\\H < ||no||«e(^(^-^-["°])-^)*. 

For r > p, we have Mp[no] < C||no||-^, so for ||no|| small enough, we have f{Mp[uo]) < fj, and the 
exponential convergence occurs. 
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When /(O) > /i, we have seen in the proof of Theorem 13.11 that L{W, Z) — )• +00 when W or Z tends 
to zero. So the trivial steady state is unstable. 



Stability of nontrivial steady states. 

Let (Woo, Zoo) be a positive steady state to System (f39l) . We want to prove that ZooV({fi; •) is locally 
asymptotically stable. Since Theorem 13.11 ensures the convergence of any solution to Equation ([8|) 
toward a steady state, it only remains to prove the local stability of ZooU{fi; •), namely 

Vpi>0, 3p2>0, \\uo- ZooU{fi;-)\\ < p2 ^ > 0, \\u{t, ■) - ZooUifi; ■)\\ < pi. 

We have already seen that 

\\u{t, •) - ZMf^; Oil < C\\qq\o - U{fi; oil + \\Q{t)ll{W{t)fi; •) - Z^U{fi; ■)\\ 

with WQUiWfi;-) - ZooUip;-)\\ ^ when {W,Z) ^ (l^oo,^oo). 
Let first treat the term ||i?Q "^uq — Oil- We have 



\qo- Zo 



{uq - ZooU {p; x))(p{fi; x) dx 

< i-r'' (^J{uo- ZooU{p;x)f{x + x'')dx 
= C\\uo- ZooK{fi--)\\ 



X 



X + x^ 



■ dx 



(49) 



and then 



Qo^\Qo - Z^\\\U{p;-)\\ 



C 



< —\\uo-ZooU{p;' 



Qo 



< 



< 



C 



Zoo ~ l^'O ~ Zoo\ 

C\\uq - ZooK{fi; 



Uq - ZooU (^; 



Zo 



C\\uq - ZooU{n\ 



(50) 



so ||^q"^wo — ^(a*; Oil is small for \\uq — ZooU{ijl\ ■)|| small enough. 

Now let turn to the term WQUiW ■) — Z^ol^ip', •)\\- Since it tends to zero when {W,Z) tends to 
(1Voo,-^oo), it is sufficient to prove that {Woo, Zoo) is stable for System (ji6]) . For t] > 0, denote 
by Vr^{Woo,Zoo) the connex component of {(W,Z), L{W,Z) < L(Woo,Zoo) + ??} which contains 
(Woo, Zoo)- Since /'(/oo) < 0, (Woo, Zoo) is a strict local minimum of L so we have, denoting B{X, p) = 

{Y E m2, ||x -y|| < p}, 



and reciprocally 



Vp>0, 3 7?>0, V^(T^oo,^oo) C5((T^oo,^oo),p) 



Vr/>0, 3p>0, B{{Woo,Zoo),p) CVr,(Woo,Zo 



So it is sufficient to have the local stability of (Woo, ^00) to prove that V,;(Woo,^oo) is stable. This 
is true for System ()39p since in this case L is a Lyapunov functional. Then, by continuity of /, there 
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exists En such that VrjiWoo, Zoo) remains stable for System ([16]) if |e(t)| < for all t > 0. But we 
know from (|45]) that \e{t)\ < C\\qq^uo —U{fi;-)\\ and from (f50]) that ||^)Q^no —U{iJ.;-)\\ is small for 
||uo — ZooU{^; ■)\\ small enough. Finally \\QUiW fi; ■) — Z^cl^ifJ-', Oil small for \\uo — ZoJJ{^; •)|| small 
enough and the local asymptotical stability of the nontrivial steady states is proved. 

Now assume that k = 2 and prove the local exponential stability of ZooU{fi; •). In this case we have 
(see examples in [23]) the explicit formula 

U{x) = Ce~^'''' (51) 
and thanks to this we can estimate the quantity \\Ql/({Wfi; •) — ZooU{fj,; •)||. We have 

WQUiW^;-)- ZooU{ti;-)f = j \QU{W'' x) - Z^U{fi; x)Wx + x') dx 

< \Z - Z^\^W-^''P Ju{W''fi;xf{x + x'-)dx (z) 
+|iy-'^(p+i) - l|2^^ J U{W-''n-''x)\x + x'-) dx {ii) 

+Z^ j \U{W-^^l-^x) -U{^l-^x)\^{x + x'')dx (m) 

and we prove exponential decay of (z), (ii) and (iii) in a neighbourhood of (Woo, Zoo)- We have thanks 
to the L'Hopital rule 

Z^Zoo —2Zoof' (Zoo) 

and 

W^Woc 2 Woo 

SO the following local "entropy - entropy dissipation" inequality holds 

3C>0,p>0, y{W,Z) e B{{Woo,Zoo),p), L{W,Z)<CD{W,Z). (54) 

Fix such a p and fix > such that VrjiWoo, Zoo) C B {{Woo, Zoo, p)- Consider ||tio — Zool^{p',-)\\ 
small enough so that |e(t)| remains smaller than for all time. Then Vr^{Woo, Zoo) is stable for the 
dynamics of System dM]). Now look at the term E{W{t), Z{t),e{t)) for (W, Z) solution to System (06]) 
in this stable neighbourhood. It satisfies 

E = p{W-l){U{{l+E)Z)-fp{Z)) + {al+al)p{fi-fj,{Z)){fj,{{l + e)Z)-fp{Z)) 

< ^AW - I)' + ^alp{p - U{Z))f + (^^^)^) (/^((l + ^)^) - fviZ))' 
= ^D{W,Z) + C{U{{l + e)Z)-fp{Z))^ 

< lD{W,Z) + Csnp\f'\e^ 
2 J 

where J = [Zoo — P — ^ri. Zoo + P + £■/?]• As a consequence we have 

^L(iy,Z) < -}-D{W,Z)^Ce^ 

< -CL{W, Z) + C\\q^\o - U{p- .)||2e-2«Mt) 
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and so 

L{W, Z) <L{Wo, Zo)e-^* + C7|| - U{^l■ •) fe'^-Mt) . 

Then, thanks to the equivalences (j52p and ()53p and because h[t) ~ t when t — )• oo, there are 6 > 
and C > such that 

(T^ - 1)2 + - < C{^l - U{Zo))^e-''" + CWq^'uo - ^(/^; ■)fe-''' 

and, because /'(/oo) / 0, it holds 

{W - 1)2 + (Z - Zoo)2 < C{Zo - Z^fe~^'' + C\\uo - Z^U{fi; Of e'^^*. 

Using the Cauchy-Schwarz inequality as in ()49p . we find that 

{Zo-Z^f <C\\uo-Z^l({fi;-)f, 



so 

{W - 1)2 + (Z - Zoo)2 < C||no - Zoo^/(/i; ■)fe-^'K 

This inequality ensures that the terms (i) and (ii) decrease to zero exponentially fast. It only remains 
to prove the same result for {in), and for this we use the explicit formula ([5T]l . We obtain 

2 



[ ( -^(w-k^)-, _M^7 _jL(i+H/-i)a:7\2 

= C / I e ' + e T'' — 2e T'^ ^ ^ I (x + x j ax 

/2/3 -Y , , , f 2B 

e'^y {W^y + W^y'')W^dy + / e"^"^ (x + x'') dx 



where 



2/3 ^ (' 2/3 

Ci:i{W) / e"-'' x(ix + C^,.(W^) / e"^"^' x^ dx 



V;,(Vr) := + 1-2 ) 

Thanks to a Taylor expansion, we find that, locally, 

\Mw)\<c{w -if 

and so 

/>oo 

V^^'^;) - ^(^"''x))^(x + x"^) dx < C\\uo - Z^Uifi; ■)fe~'^''\ 







Finally there exists a consntant C > such that, for ||mo — ZooU{fi; ■)\\ small enough so that {W,Z) 
stay in the neighbourhood VrjiWoo, Zoo), we have 

\\u{t,-)-ZMfi;-)\\ < \\u{t,-)-QU{Wfi;-)\\ + \\QU{Wfi;-)-ZMfi;-)\\ 
< C\\uo-ZMfi;-)\\e-'" 

and it is the local exponential stability of the nontrivial steady states which satisfy /'(/oo) < in the 
case when k = 2. 

When /'(/oo) > 0, the steady state (Woo, ■^oo) is a saddle point of L so it is instable. □ 
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We remark that the structure of the reduced system (j39p is different for p < 1 and p > 1. The 
nontrivial steady states are focuses in the case when p < 1 and nodes for p > 1 (see Figure [T] for a 
numerical illustration in the case of Corollary I1.3P . 



N 




Figure 1: Solutions to System (j39p are plotted in the phase plan {W,Z) for two different values of 
parameter p. The other coefficients are 7 = 0.1, /x = 1 and /(x) = 2e~^. We can see that the steady 
state is a focus for p < 1 (left) and a node for p > 1 (right). 



4 Nonlinear Drift and Death Terms: Stable Persistent Oscillations 

We have seen in Theorem 11.11 that any solution to the nonlinear equation ([8]) converges to a steady 
state. Can this result be extended to Equation (llOp where the death rate is also nonlinear ? The 
result in Theorem 11.41 answer negatively to this question. Indeed it ensures the existence of functions 
/ and g, and parameters p and q such that Equation (jlOp admits periodic solutions. Here we give 
examples of such functions and parameters and more precisely we prove, thanks to the Poincare- 
Bendixon theorem, that any solution with an initial distribution in the eigenmanifold £ which is not 
a steady state converges to a periodic solution. Then we extend this result by surrounding this set of 
initial distributions by an open neighbourhood in T-L. 

In the proof, we need to know the dependency of some quantities on the parameters p and q. Since 
we do not know the dependencies of Mp = J x'pU{x) dx on p, we consider an equation slightly different 
from (llOp . namely 

Clearly the existence of functions / and g for which persistent oscillations appear in Equation (j55p 
ensures the same result for Equation (|10p (up to a dilation of / and g). Now let set the assumptions 
on the two function / and g which allow to obtain periodic oscillations. We consider differentiable 
increasing functions / and g such that 

/(O) > £?(0) = and /(oo) < 5(00) = 00. (56) 
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Moreover assume that g is invertible and define on M"*" the function 

^(W) := f (w''^P-'''>g~\W)'^ . (57) 
To ensure the existence and uniqueness of a nontrivial equilibrium, we assume that 

3!l^oo > 0, i'iWoo) = Woo and moreover ^p'iWoo) < 1- (58) 
This steady state is unstable if, denoting Q^o '■= WoJ"^ (Wcx:,) , we have 

Qoo (pW^PfiW^PQoo) - W^^g'iW^^Qoo)) - ^ > 0. (59) 

Then sohitions to Equation (j^^ with initial distribution close to the set £ \ {Qoo^(^oo; ")} exhibit 
asymptotically periodic behaviors. More precisely, we have the following result 

Theorem 4.1. Consider increasing dijjerentiable functions f and g satisfying conditions (j56p . (j58|) 
and ([59]) . a parameter 7 G (0, 2] and a fragmentation kernel k, which satisfies Assumption ([32]) . Then 
there exists an open neighbourhood V of £\{Qao^{Woo', •)} ^ such that, for any initial distribution 
uq G V, there exists periodic functions W{t) and Q{t) such that 

\\u{t, •) - Qit)U{W{ty, Ollw y 0. (60) 

Before proving this theorem, we give examples of functions / and g such that conditions (j56p . (j58|) 
and (I59p are satisfied. These examples, together with Theorem 14.11 give the proof of Theorem 11.41 

Example 1. Assume that there exists C > such that for all a; > 0, g{x) < C xg'(x), then ^lJ(W) = W 
has a unique solution for k(q — p) > C. Indeed, if we compute the derivative of we find 

^'(VT) = T^Mp-.) (^k^p _ g)^M + (g-^Yiw)^ 

and g{x) < C xg'{x) implies that x{g~^y{x) < C g~^{x). So if k{q — p) > C, tp decreases and Assump- 
tion ()58p is fulfilled. If moreover /(I) = g{l) = 1, then the unique nontrivial equilibrium is given by 

Woo = 1- Then condition ([59]) is satisfied for p > ^-j^^^. 

Example 2. Consider the case g{x) = x and p = q, and assume that f{x) — x has a unique root xq 
and /'(xq) — 1 < 0. Then il){W) = fiW) and Assumption ([^5j) is satisfied. Moreover, condition ([^^ 
writes pf'{Woo) > 1 + ^ so it is satisfied for p large enough. 

Now we give a lemma useful for the proof of Theorem 14.11 
Lemma 4.2. Consider a dynamical system in M" with a parameter e{t) 

X = F{X-e{t)) (61) 

with F € C(M" X M"). Assume that for any vanishing parameter ||e(t)|| *~^°^> the solutions to 
Equation (j6ip are hounded. Then for any solution associated to ||e(t)|| *~^°°> 0, there exists a 
solution X^ associated to e = such that X^ and X^ have the same lo— limit set. 
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Proof of Lemma \4 ■ 4 Let X{t) a solution to System ([6T|) with ||e(t)|| — )■ 0. By assumption, X{t) is 
bounded, so X{t) is also bounded since F is continuous. Now consider a sequence {tfcjfceN wich tends 
to the infinity and define the sequence {VFfc(-)} by Wk{t) = W{t + t^). This sequence is bounded 
in VF-'^'°°(M_|_) so there exists a subsequence which converge to Woo{-)- This limit is a solution to 
Equation (f6T]) with e = 0. We take := Woo that ends the proof of Lemma HT2l □ 



Proof of Theorem \4-l\ We divide the proof in two parts: first the result for mq G <f and then the 
existence of a neighbourhood V of £ in T-L where the result persists. 

First step: uq G £. 

For uq £ £ \ {0}, there are Wq > and Qq > such that 

uoix) = Qol^iWo;x). 
Then, if u(t, x) is the solution to Equation ()55p and W is the solution to 

W ( ^ (Mp[u]{t) 



with W(0) = Wo, the relation holds for ah t > and x > 

u{t,x) = Q{t)U{W{t)-x) 
where Q{t) := QQe-l'o(^i^)-9i^M^]i^)/M,[mds ^ Tj^^^^ compute 

/■oo 

MpH(t) = / xPu{t,x) dx = W^P{t)Q{t) Mp[U] 







and finally we obtain the reduced system of ODEs satisfied by {W, Q) 

W . . _ , 

W 



' W = — (/ {w^PQ 

Q = Q(w-g(w'"^Q)). 



(62) 



We prove that System (j62p has bounded solutions and a unique positive steady state which is unstable. 
Then we use the Poincare-Bendixon theorem to ensure the convergence to a limit cycle. 

The fact that < /(O) < / < /(oo) < oo and that g increases from to the oo ensures that the 
solution remains bounded. Let {Woo,Qoo) a positive steady state. It satisfies 



Woc = f (W'pq) = g (w'^^q' 
and so, since g is invertible, Qoo = W^'"^ g~^ (Woo) ■ Then Woo is solution to the equation 

Woo=f (^4(^-^)5-1 (Woo)) = V'(W^oo) 
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but Assumption (|58p ensures the uniqueness of such a solution. Now look at the stability of this 
positive steady state. We write system (I62p under the form 

W\ (W 

Q ) \Q 

and we have 

J (F\ =( vW^Q^f{W^Q^)-lW^ IW^^'nW^Q^) 
The trace of this matrix is 



k 



and the determinant is 



D = (w'"ig\W'"^Q) - W''Pf'{W''PQ)^ + {q - p)W''^P+i^Q^f'{W''PQ)g'{W'"'Q). 
We know from Assumption (j58|) that ip'{Woo) < 1 and, if we compute ip'{W) we find 



Since g '^{W) = W^'^Q we finally obtain 



D = \wt^^^Q^ g'{W^^Q^){l - ^'(t^oo)) > 0. 



Thus when T > 0, namely when Assumption (159p is satisfied, the two eigenvalues have positive real 
parts and the positive steady state is unstable. Now we prove that iyV^ Q) remains away from the 
boundaries of (M+)^. For this we write that 

Vt > 0, W_ := min(VFo, /(O)) < W{t) < max(Wo, /(oo)) := W 

and then 

Q>mm{Qo,W~''%~\W)). 

Since f{0)>0,W_>0 for Wq > so any solution with Wq > and Qo > stays at positive distance 
of the boundaries of (M_|_)^. Then the Poincare-Bendixon theorem (see [35] for instance) ensures any 
solution to System (l62]) with Wq > 0, Qq > and {Wq, Qq) ^ {Woo, Qoo) converges to a limit cycle. 

Second step: Existence of V. 

Let uq ^ in Ti and define from u{t, x) solution to Equation (j55p with initial distribution uq a function 
V by 

v{h{t),x) = VF'=(t)n(t,VF^Xi)a;)e^o^'^^*^^f"'^'^/*^''"^^"^^'^^'^' 

with W solution to 
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and h solution to h = W with h{0) = 0. We have already seen in Section [2^31 that h : M+ — is one 
to one since h{t) > kln{l + |). We take W{0) = 1 to have v{t = 0, •) = u{t = 0, •) = uq. Thanks to 
Theorem 12.11 we know that v is solution to 

dtv{t, x) + dx{x v{t, x)) + v{t, x) = T^v{t, x) 

and the GRE ensures the convergence 

v{t,x) > I / (I){x)uq{x) dx ] U{x). 

As a consequence we have the equivalences, for any p > 0, 

Mp[u]{t) ~ gf^Mp[U]W''Pe^oiW{s)-g{M,M{s)/M,[U])ds 

so, if we define Q{t) := ^oe-/'o(^(*)-5(*^9N(«)/A'/,[W])rfs^ that the reduced system ([62]) is "asymp- 

totically equivalent" to Equation (|55p . More precisely we define as in the proof of Theorem 13. II 

Spit) - 



Mp[U]Wf'P{t)Q{t) 
we have that — t- and {W, Q) is solution to 



Q = Q(W-gi{l + e, 



-1 



(63) 



Now we prove that if Wq and Qq are positive and (Wq) Qo) / (W^ooi Qoo)) then if ||tio — QqU{Wq; •)||-^ 
is small enough, the solution u to Equation ([55]) converges to a periodic solution. Denote by d the 
distance between (Wo,Qo) ^-iid (WqcQoo)- Since (Woo,Qoo) is a source for System there exists 
a ball with radius p < d such that the flux is outgoing, namely 

y{W,Q)€dB{{W^,Q^),p), F{W,Q)- n>0 (64) 

where n is the outgoing normal of B{{Woc,Qoc), p)- Then if we define by F{W,Q;ep,eg) the flux of 
Equation (j63p . we have by continuity of / and g that there exists Eq such that (|64p remains true for 
-F(VF, Q; Ep, Eg) provided that and £q stay less than eg. But we know from the proof of Theorem 13. II 
that there exists a constant Cp > such that for all time t > 0, Ep < Cp\\uo — QoU{Wo; •)\\- So for 
ll^^o — Qo^iWo; •)|| < Q-^f^, the solution to System ([63]) cannot converge to the positive steady state 
(Woo, Qoo)- Thanks to the same arguments, if \\uo — QoU{Wo; •)|| is small enough, then (W, Q) remains 
away from the boundaries of (M+)^. We obtain thanks to Lemma [4 . 2 1 that for ||no — QqIA{Wq; •)|| small 
enough, iW {t),Q{t)) converges to a limit cycle {W {t),Q{t)). Then we write 

\\u - QUiW; Oil < ||n - QKiW; ■)\\ + ||Q^(P^; •) - QU{W; ■)\\ 

and we conclude as in the proof of Theorem 13.11 that the solution u to Equation \W\ converges in to 
QU{W; ■). Finally we have proved, for any {Wq,Qq) £ (R^)'^ \ {(Woo, Qoo)}, the existence of a ball 
centered in QqU{Wo; •) such that any solution to Equation (|55p with an initial distribution in this ball 
converges to a periodic solution. Then Theorem 14.11 is proved for V the reunion of all these balls. 

□ 



23 



To illustrate the convergence to a periodic solution for solutions to Equation (j55p . we plot in Figure [2] 
a solution to Equation ([62]) with an initial condition close to the steady state {Woo,Qoo) and for 
coefficients which satisfy the assumptions of Theorem 14.11 




5 The Prion Equation: Existence of Periodic Solutions 



Prion diseases are believed to be due to self-replication of a pathogenic protein through a polymer- 
ization process not yet very well understood (see [40] for more details). To investigate the replication 
process of this protein, a mathematical PDE model was introduced by [32]. We recall this model here 
under a form slightly different from the original one (see [121 [20] for the motivations to consider this 
form) 

dV{t) 



dt 



X - V{t) 6+1 T{x)u{t,x) dx 
I Jo 

d 

—V(t) — (T(x)u(t, x)) — fi(x)u(t, x) + Fuit, x) 
ox 



(65) 



In this equation, u(t, x) represents the quantity of polymers of pathogenic proteins of size x at time t, 
and V{t) the quantity of normal proteins (also called monomers). The polymers lengthen by attaching 
monomers with the rate r(x), die with the rate ijl{x) and split into smaller polymers with respect to 
the fragmentation operator The quantity of monomers is driven by an ODE with a death parameter 
5 and production rate A. This ODE is quadratically coupled to the growth-fragmentation equation 
because of the polymerization mechanism which is assumed to follow the mass action law. 

This system admits a trivial steady state, also called disease-free equilibrium since it corresponds 
to a situation where no pathogenic polymers are present: V = ^ and u = 0. The stability of this 
steady state has been investigated under general assumptions on the coefficients in (TTJ [T21 ESI ES] ■ 
The results indicate that the disease-free equilibrium is stable when it is the only steady state, and 
becomes unstable when a nontrivial steady state appears (also called endmic equilibrium). It is proved 
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in [9] that several such nontrivial steady states can exist. But the stabihty (even hnear) of these steady 
states is a difficult and still open problem for general coefficients. The only existing results concern 
the "constant case" (r constant, (3 linear and k constant) initially considered by [32], since then the 
model reduces to a closed system of ODEs. In this case, the problem has been entirely solved by 
[241 [32| I55j: the endemic equilibrium is unique and it is globally stable when it exists. 

A new more general model has been introduced in [S^ and takes into account the incidence of 
the total mass of polymers P{t) := J xu{t,x) dx on the polymerization process. More precisely, they 
consider that the presence of many polymers reduces the attaching process of monomers to polymers by 
multiplying the polymerization rate by x^^Jp^f) with uj a positive parameter. Then they prove similar 
results about the existence and stability of steady states, still in the case of constant parameters. 

Here we look at a generalization of the influence of polymers on the polymerization rate by consid- 
ering the system 

^^Ki^l = -V(t)f ( [ xPu] r T(x)u(t,x) dx-6V(t) + X, 

)^ (66) 
—u{t,x) = -V{t)f i xPuj —{T{x)u{t,x)) - fiu{t,x) +Tu{t,x), 

where p > and / : — )• M"*" is a differentiable function. In this framework, the model of [33] 
corresponds to p = 1 and f{P) = i_^}^p , together with r constant, /3 linear and k constant. Using the 
reduction method to ODEs, we prove that such a system can exhibit periodic solutions. For this we 
consider the foolowing system, where Mp = Mp^], 

'^^^^^ - -V{t)f (fi^'^P^ [ x^u] r X u{t, x) dx - 6V{t) + A , 



dt V M, 



p J y Jo 



(67) 

—u{t, x) = -V{t)f J x^uj — {x u{t, x)) - fi u{t, x) + J^^uit, x), 

which is a particular case of System (|66p with coefficients satisfying the assumptions of Theorem 12.11 
and up to a dilation of /. We prove that, under some conditions on the incidence function /, there 
exist values of p for which System (j67p admits periodic solutions. The method is to consider p as a 
varying bifurcation parameter and prove that Hopf bifurcation occurs when p increases. 

Theorem 5.1. Define on [0, ^) the function 



A — /i'^+^x 

and consider a positive differentiable function f satisfying 

3\xo > s.t. f{xo) = g{xo) and moreover < f' {xq) < g{xQ). (68) 

Assume that fi < {k + fi~^)5, then there exists p > for which Equation (|67|) admits a periodic solution. 
More precisely, there exist V, W and Q periodic such that 

{V{t),Q{t)U{W{t);x)) is solution to Equation ([67]). 



25 



Proof. First step: reduced dynamic in £ 

We look at the dynamic of System (167p on the invariant eigenmanifold £. For any initial condition uq 
in £, there exist Qq and Wq such that uq writes 

uoix) = ■^QoKiWo;x). 

Consider the solution to System ([67|) corresponding to this initial data and define W as the solution 
to 



W 



k 



W =—{f{ ^-^P^ I xPu]V-W 



Mr, 



(69) 



_ W{^) = Wo. 
Then we know from Theorem 12.11 that the solution satisfies 



and this allows to compute 



and 



u{t,x) = 9^U{W{t)-x)e^oms)-^^)ds 



/ xPu(t, x) dx = Qo^W^Pe^o iw{s)-^l) ds 
Jo Ml 



X u{t, x) dx = QqW^ eio iW[s)-f^) ds_ 
Thus, if define Q{t) := Qo e-^o Equation §9i) becomes 



W = ^(^f (^{fi^^W)^PQ )V-W 



and System (|67|) reduces to 



V = x-v (5 + f ({h^^w)^pq)w''q 



W 



(70) 



I Q = Q{w-f,). 



Now we prove that System (jTOp admits a unique nontrivial steady state which undergoes a supercritical 
Hopf bifurcation when p increases from 0. 

Second step: Hopf bifurcation for the reduced system 

First we look for a positive steady state of System (j70|) . Such a steady state is unique and given by 



1 



A - /+iQ, 



where Qoo satisfies 



fiQc 



6fj, 



=: giQc 
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Such a Qoo exists and is unique by Assumption ([68]) . moreover it satisfies A — fi^^^Qoo > and so Voo 
is positive. Finally there exists a unique positive steady state and we prove that a Hopf bifurcation 
occurs at this point. The linear stability of the steady state is given by the eigenvalues of the Jacobian 
matrix 

-fc/FQ(/(Q)+pQ/'(Q)) -/F(/(Q) + Q/'(Q)) 
Jace, = ^f{Q) Pf^VQf'iQ) - f IVf'iQ) 

V Q 

where the oo indices are suppressed for the sake of clarity. The trace of this matrix is 

T = -6- fi'^QfiQ) - ^ + PisVQfiQ) 

which is negative for p < pi and positive for p > pi with 

^ S + fi/k + fi'^QfiQ) 

The determinant is 
It is independent of p and negative since /'(xq) < g'{xQ) and 

(A — fi''+'-xo) 
The sum of the three 2x2 principal minors is 

M = -SpVQf'iQ) + ^ + /(/z + ^)Q/(Q) - ^VQf'iQ). 

To use the Routh-Hurwitz criterion, let define ip{p) := MT — D and look at its sign. For p = we 
have 



+ 5Qf{Q) + ^+ {5{k + ^-1) - /i) /g/(Q) 



+ VQ (/-i(A; + - nQ)) (l^'QfiQ) + T 



and it is negative since ^ < (k + ^x-^)8 and f'{Q) < g'{Q) = ^f{Q) < fi''-^{k + fi-^)f{Q). For 
p = pi, it is positive because ip{pi) = —D > 0. Now we investigate the variations of ip between and 
Pi. The first derivative of T, D and M are given by 

T'{p)=t,VQf'{Q), M'{p) = -6VQf'iQ), D'{p) = 0, 

and the second derivatives are all null 

T"{p) = M"{p) = D"{p) = 0. 

So we have 

= 2M'{p)T'{p) < 
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and is concave. Thus there exists a unique po € (0,pi) such that ip{po) = 0. Now we can use the 
Routh-Hurwitz criterion (see |35j for instance). For < p < po we have T < 0, D < and MT < D, 
so the steady state is linearly stable with one real negative eigenvalue and two complex conjugate 
eigenvalues with a negative real part. For pQ < p < pi we have T < 0, D < and MT > D, so 
the steady state is linearly unstable with one real negative eigenvalue and two complex conjugate 
eigenvalues with a positive real part. The two conjugate eigenvalues cross the imaginary axis when 
p = Po so there is a Hopf bifurcation at this point. To prove that a periodic solution appears with this 
bifurcation, it remains to check that the complex eigenvalues cross the imaginary axis with a positive 
speed (see [29] for instance). Denote by a it z6 the two conjugate eigenvalues and c < the real one. 
We have to prove that the derivative a'{po) > 0. For this we express ip{p) in terms of a{p), h{p) and 
c(p), and we use the concavity of Tp. We have for any p 

T = 2a + c, D = c{a^ + b'^), M = a^ + b^ + 2ac, 

so 

ipip) = 2a{a^ + h^) + Aa^c + 2ac^. 
Then, using that a{po) = by definition of po, we obtain 

i;'{po) = 2{b^ + c^)a'. 

But Tp'ipo) > because ■0 is concave and increasing on a neighbourhood of po, so necessarily a'(po) > 0. 
This proves the existence of a periodic solution {V, W, Q) to System ([70]l for a paramter p > po close 
to pq. Then the functions V{t) and Q{t)U{W{t),x) are periodic and solution to System ([67j) . □ 

The question to know if such a periodic solution is stable is difficult, even for the reduced dy- 
namic (j70p . Nevertheless we give in Figure [3] evidences that it should be the case. This simulation is 
made with parameters and a function / satisfying Assumption (I68p . for a value of parameter p > pQ. 
It seems to indicate that the periodic solution persists for p away from po- 
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6 Comparison between Perron and Floquet Eigenvalues 



In this Section, we consider that the time-dependent terms V{t) and R{t) of the growth-fragmentation 
equation are T-periodic controls. Periodic controls are usually used in structured equations to model 
optimization problems. In the case of prion diseases (see Section [5]), there exists an amplification 
protocole called PMCA (see [IQ] and references therein) which consists in periodically sonicating a 
sample of prion polymers in order the break them into smaller ones and thus increase their quantity. 
Between these phases of sonication, the sample is let in presence of a large quantity of monomers 
in order to allow a fast polymerization process. This protocole can be modeled by introducing in 
the growth-fragmentation equation a periodic control in front of the fragmentation operator [9l llOj . 
Then a problem is to find a periodic control which maximizes the proliferation rate of the polymers 
in the sample. Mathematically it can result into the problem to optimize the Floquet eigenvalue of 
the growth-fragmentation equation, namely the eigenvalue associated to periodic coefficients. Before 
solving this difficult question, a first step is to compare the Floquet eigenvalue to the Perron eigenvalue 
associated to constant coefficients, for instance the mean value on a period of the periodic ones, and 
to know if the Floquet one can be better than the Perron one. Such concerns are also investigated 
in the context of circadian rhythm for the optimization of chronotherapy (see \13\ HM I15|). The 
population is an age structured population of cells and the model is a system of renewal equations. 
The death and birth rates are assumed to be periodic, and the Floquet eigenvalue is compared to the 
Perron eigenvalue associated to geometrical or arithmetical time average of the periodic coefficients. 
Comparison results obtained show that the Floquet eigenvalue can be greater or less than the Perron 
one depending on parameters. 

Here the controls are on the growth and death coefficients, and we give comparison results between 
Floquet and Perron eigenvalues in the case where ;y = 1 or z/ = and 7 = 1. The Floquet eigenelements 
{KfMf) associated to periodic controls are defined by two properties: UF{t,x)e^^^ is a solution to 
Equation ([12]) and Upit^x) is a T-periodic function of the time. For any T-periodic function /(t), we 
use the notation 



To ensure the uniqueness of Floquet eigenfunction, we imposes Up = ^- Then we have the following 
comparison results. 

Proposition 6.1. Assume that conditions ©-([T]) are satisfied with v = 1, and define, from V{t) the 
T-periodic control on growth, W{t) as the periodic solution to 




W = 



tW 



(V-W). 



(71) 



k 



Then the identities hold 




(72) 



and 



Af = A{V,R) = A{V,R). 



(73) 



Proof. Thanks to Corollarv 12.31 we have that 



U{Wit);x)e^oHW{s)Ms))ds 



Upit, x)e 



Apt 
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oo 



and so, integrating in x, we find 
Then, by periodicity, we have 







^ApT-f^ A(W{s)Ms))ds ^ ^ 

which gives 

r-T 1 rT 



^F = ^l^ A{W{s),R{s))ds = ^ 1^ {TW{s)-fiR{s))ds. 



Thanks to ODE ([Til) we have 

rT 

V-W = 



and so 

rT 







Af = ^^ {TV{s)-iiR{s))ds = ^ A{V{s),R{s))ds. 



□ 



In the case i/ = and 7 = 1, we cannot ensure the existence of Floquet eigenelements with our 
method. Nevertheless, thanks to CoroUary 12.51 we can compare the eigenvalues of the reduced sys- 
tem (j26p which is satified by Mo[ii] and Mi[n]. 

Proposition 6.2. Assume that t{x) = r, /3(x) = (3x, fi{x) = ji and n is symmetric, then we have the 
comparison 

A{V,R) <Af< A{V,R). (74) 
Proof. Define W as the periodic solution to 

w = Mi^iv - W). 

k 

We know thanks to CoroUarv 12.51 that 

Mo[W]{t) = Mo[U]e^oMW{s),R{s))ds 

and 

Mi[W]{t) = Mi\U]W''{t)efoMW{s),R{s))ds 
are solution to System (p6|) . As a consequence 

^F = ^£ AiWis), R{s)) ds = i j\^pTW{s) - fiRis)) ds. 
Using the ODE satisfied by W, we have 

F V -W 
'0 W~ ~k~ Jo ~^/W 
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and we obtain that 

rT y 



W 







W 



Then the Caucy-Schwartz inequahty gives 



Jo W^ \Jo Vw\Jo ' 



and so 

fT ^ rT i rT 



To obtain the second inequahty in (I74p . we write using the ODE satisfied by W 



w vtt r 

0= / = ^4— / V-W 



and so 



Thus we have, using the Jensen inequahty, 



w 




Vw 


k 








-f 


Jo 


Jo 



^ ' yJW{s)ds < [ W{s)ds = J^f V{s)ds 



and finaUy 



Af = ^^ iVPrW{s)-fiR{s))ds<J^j^ PTV{s)ds-^j^ lJiR{s) ds = K{V ,R) . 



□ 



Conclusion and Perspectives 

We have introduced a new reduction method to investigate the long-time asymptotic behavior for 
growth-fragmentation equations with a nonhnear growth term. It ahows to prove convergence and 
stabihty results but also to prove the existence of periodic solutions by using ODE methods. The 
technique is based on self-similar dependencies which require powerlaw coefficients, and the main 
theorem requires moreover a linear size-dependency of the growth coefficient. Nevertheless, we have 
given some results about the case when the growth rate is only assumed to be a powerlaw. A step 
further would be to adapt the ideas in this paper to build nonlinear entropy techniques adapted to 
general powerlaw growth coefficients. 



Aknowledgment 

The author would like to thank Jean-Pierre Erangoise for helpful discussions about limit cycles. 



31 



References 

[1] M. Adimy, F. Crauste, M. L. Hbid, and R. Qesmi. Stability and hopf bifurcation for a cell popula- 
tion model with state-dependent delay. SIAM Journal on Applied Mathematics, 70(5):1611-1633, 
2010. 

[2] H. Banks, F. Charles, M. Doumic, K. Sutton, and W. Thompson. Label structured cell prolifer- 
ation models. App. Math. Letters, 23(12):1412-1415, 2010. 

[3] V. Bansaye and V. C. Tran. Branching feller diffusion for cell division with parasite infection. 
ALE A, to appear. 

[4] F. Bekkal Brikci, J. Clairambault, and B. Perthame. Analysis of a molecular structured population 
model with possible polynomial growth for the cell division cycle. Math. Comput. Modelling, 47(7- 
8):699-713, 2008. 

[5] F. Bekkal Brikci, J. Clairambault, B. Ribba, and B. Perthame. An age-and-cyclin-structured cell 
population model for healthy and tumoral tissues. J. Math. Biol., 57(1):91-110, 2008. 

[6] S. Bertoni. Periodic solutions for non-linear equations of structured populations. J. Math. Anal. 
Appl, 220(1) :250-267, 1998. 

[7] T. Biben, J.-C. Geminard, and F. Melo. Dynamics of Bio-Polymeric Brushes Growing from a 
Cellular Membrane: Tentative Modelling of the Actin Turnover within an Adhesion Unit; the 
Podosome. J. Biol. Phys., 31:87-120, 2005. 

[8] M. J. Caceres, J. A. Cahizo, and S. Mischler. Rate of convergence to an asymptotic profile for 
the self-similar fragmentation and growth-fragmentation equations. To appear in J. Math. Pures 
Appl, 2011. 

[9] V. Calvez, M. Doumic Jauffret, and P. Gabriel. Self-similarity in a general aggregation- 
fragmentation problem ; application to fitness analysis, preprint, arXiv:1011.5468, November 
2010. 

[10] V. Calvez and P. Gabriel. Optimal control for a discrete aggregation-fragmentation model, in 
preparation. 

[11] V. Calvez, N. Lenuzza, M. Doumic, J. -P. Deslys, F. Mouthon, and B. Perthame. Prion dynamic 
with size dependency - strain phenomena. J. of Biol. Dyn., 4(l):28-42, 2010. 

[12] V. Calvez, N. Lenuzza, D. Oelz, J. -P. Deslys, P. Laurent, F. Mouthon, and B. Perthame. Size 
distribution dependence of prion aggregates infectivity. Math. Biosci., 1:88-99, 2009. 

[13] J. Clairambault, S. Gaubert, and T. Lepoutre. Comparison of Perron and Floquet eigenvalues in 
age structured cell division cycle models. Math. Model. Nat. Phenom., 4(3): 183-209, 2009. 

[14] J. Clairambault, S. Gaubert, and T. Lepoutre. Circadian rhythm and cell population growth. 
Math. Comput. Modelling, doi:10.1016/j.mcm.2010.05.034. In Press, 2010. 

[15] J. Clairambault, S. Gaubert, and B. Perthame. An inequality for the Perron and Floquet eigen- 
values of monotone differential systems and age structured equations. C. R. Math. Acad. Sci. 
Paris, 345(10):549-554, 2007. 



32 



[16] J. M. Gushing. Bifurcation of time periodic solutions of the McKendrick equations with apph- 
cations to population dynamics. Comput. Math. AppL, 9(3):459-478, 1983. Hyperbolic partial 
differential equations. 

[17] O. Destaing, F. Saltel, J.-C. Geminard, P. Jurdic, and F. Bard. Podosomes Display Actin 
Turnover and Dynamic Self- Organization in Osteoclasts Expressing Actin-Green Fluorescent 
Protein. Mol. Biol, of the Cell, 14:407-416, 2003. 

[18] O. Diekmann, M. Gyllenberg, and J. A. J. Metz. Steady-state analysis of structured population 
models. Theoretical Population Biology, 63:309-338, 2003. 

[19] M. Doumic. Analysis of a population model structured by the cells molecular content. Math. 
Model. Nat. Phenom., 2(3):121-152, 2007. 

[20] M. Doumic, T. Goudon, and T. Lepoutre. Scaling limit of a discrete prion dynamics model. 
Comm. Math. Sci., 7(4):839-865, 2009. 

[21] M. Doumic, P. Maia, and J. Zubelli. On the calibration of a size-structured population model 
from experimental data. Acta Biotheoretica, 2010. 

[22] M. Doumic, B. Perthame, and J. Zubelli. Numerical solution of an inverse problem in size- 
structured population dynamics. Inverse Problems, 25(electronic version) :045008, 2009. 

[23] M. Doumic Jauffret and P. Gabriel. Eigenelements of a general aggregation-fragmentation model. 
Math. Models Methods Appl. Sci., 20(5):757-783, 2010. 

[24] H. Engler, J. Priiss, and G. Webb. Analysis of a model for the dynamics of prions ii. J. Math. 
Anal. Appl, 324(1):98-117, 2006. 

[25] M. Escobedo, S. Mischler, and M. Rodriguez Ricard. On self-similarity and stationary problem for 
fragmentation and coagulation models. Ann. Inst. H. Poincare Anal. Non Lineaire, 22(1):99-125, 
2005. 

[26] J. Z. Farkas. Stability conditions for a nonlinear size-structured model. Nonlin. Anal. Real World 
Appl, 6:962-969, 2005. 

[27] J. Z. Farkas and T. Hagen. Stability and regularity results for a size-structured population model. 
J. Math. Anal Appl, 328(1):119 - 136, 2007. 

[28] J. Z. Farkas and T. Hagen. Asymptotic analysis of a size-structured cannibalism model with 
infinite dimensional environmental feedback. Comm. Pure Appl. Anal, 8:1825-1839, 2009. 

[29] J. -P. Frangoise. Oscillations en biologic, volume 46 of Mathematiques & Applications (Berlin) 
[Mathematics & Applications]. Springer- Verlag, Berlin, 2005. Analyse qualitative et modeles. 
[Qualitative analysis and models]. 

[30] P. Gabriel. The shape of the polymerization rate in the prion equation. Math. Comput. Modelling, 
doi:10.1016/j.mcm.2010.03.032. In Press, 2010. 

[31] P. Gabriel and L. M. Tine. High-order WENO scheme for polymerization-type equations. ESAIM 
Proc, 30:54-70, 2010. 



33 



[32] M. L. Greer, L. Pujo-Menjouet, and G. F. Webb. A mathematical analysis of the dynamics of 
prion proliferation. J. Theoret. Biol, 242(3):598-606, 2006. 

[33] M. L. Greer, P. van den Driessche, L. Wang, and G. F. Webb. Effects of general incidence and 
polymer joining on nucleated polymerization in a model of prion proliferation. SIAM Journal on 
Applied Mathematics, 68(1):154-170, 2007. 

[34] M. Gyllenberg and G. F. Webb. A Nonlinear Structured Population Model of Tumor Growth 
With Quiescence. J. Math. Biol, 28:671-694, 1990. 

[35] J. Hofbauer and K. Sigmund. The theory of evolution and dynamical systems, volume 7 of 
London Mathematical Society Student Texts. Cambridge University Press, Cambridge, 1988. 
Mathematical aspects of selection. Translated from the German. 

[36] F. Hoppensteadt. Mathematical theories of populations: demographics, genetics and epidemics. 
Society for Industrial and Applied Mathematics, Philadelphia, Pa., 1975. Regional Conference 
Series in Applied Mathematics. 

[37] M. lannelli. Mathematical theory of age- structured population dynamics, volume 7 of Applied 
mathematics monographs C.N.R. Giardini editor! e stampatori, Pisa, 1995. 

[38] T. Kostova and J. Li. Oscillations and stability due to juvenile competitive effects on adult 
fertility. Computers & Mathematics with Applications, 32(11):57 - 70, 1996. BIOMATH-95. 

[39] P. Laurengot and B. Perthame. Exponential decay for the growth- fragmentation/cell-division 
equation. Commun. Math. ScL, 7(2):503-510, 2009. 

[40] N. Lenuzza. Modelisation de la replication des Prions: implication de la dependance en taille des 
agregats de PrP et de I'heterogeneite des populations cellulaires. PhD thesis, Paris, 2009. 

[41] Z. Liu, P. Magal, and S. Ruan. Hopf bifurcation for non-densely defined Cauchy problems. 
Zeitschrift fur Angewandte Mathematik und Physik, doi:10.1007/s00033-010-0088-x, 2010. 

[42] P. Magal, P. Hinow, F. Le Foil, and G. F. Webb. Analysis of a model for transfer phenomena in 
biological populations. SIAM Journal on Applied Mathematics, 70:4062, 2009. 

[43] P. Magal and S. Ruan. Center manifolds for semilinear equations with non-dense domain and ap- 
plications to Hopf bifurcation in age structured models. Mem. Amer. Math. Soc, 202(951) :vi-|-71, 
2009. 

[44] P. Magal and S. Ruan. Sustained oscillations in an evolutionary epidemiological model of influenza 
A drift. Proc. R. Soc. Land. Ser. A Math. Phys. Eng. Sci., 466(2116):965-992, 2010. 

[45] J. A. J. Metz and O. Diekmann, editors. The dynamics of physiologically structured populations, 
volume 68 of Lecture Notes in Biomathematics. Springer-Verlag, Berlin, 1986. Papers from the 
colloquium held in Amsterdam, 1983. 

[46] P. Michel. Existence of a solution to the cell division eigenproblem. Math. Models Methods Appl. 
Sci., 16(7, suppl.):1125-1153, 2006. 



34 



[47] P. Michel. General Relative Entropy in a nonlinear McKendrick model. In G.-Q. Chen, E. Hsu, 
and M. Pinsky, editors, Stochastic Analysis and Partial Differential Equations, volume 429, pages 
205-232. AMS: Contemporary Mathematics, 2007. 

[48] P. Michel, S. Mischler, and B. Perthame. General entropy equations for structured population 
models and scattering. C. R. Math. Acad. Sci. Paris, 338(9) :697-702, 2004. 

[49] P. Michel, S. Mischler, and B. Perthame. General relative entropy inequality: an illustration on 
growth models. J. Math. Pures Appl. (9), 84(9):1235~1260, 2005. 

[50] S. Mischler, B. Perthame, and L. Ryzhik. Stability in a nonlinear population maturation model. 
Math. Models Methods Appl. Sci., 12(12):1751-1772, 2002. 

[51] B. Perthame. Transport equations in biology. Frontiers in Mathematics. Birkhauser Verlag, Basel, 
2007. 

[52] B. Perthame and L. Ryzhik. Exponential decay for the fragmentation or cell-division equation. 
J. Differential Equations, 210(1):155-177, 2005. 

[53] B. Perthame and S. Tumuluri. Nonlinear renewal equations. In Selected topics in cancer modeling, 
Model. Simul. Sci. Eng. Technol., pages 65-96. Birkhauser Boston, Boston, MA, 2008. 

[54] J. Priiss. On the qualitative behaviour of populations with age-specific interactions. Comput. 
Math. Appl., 9(3):327-339, 1983. Hyperbolic partial differential equations. 

[55] J. Priiss, L. Pujo-Menjouet, G. Webb, and R. Zacher. Analysis of a model for the dynamics of 
prion. Dis. Cont. Dyn. Sys. Ser. B, 6(l):225-235, 2006. 

[56] G. Simonett and C. Walker. On the solvability of a mathematical model for prion proliferation. 
J. Math. Anal. Appl, 324(l):580-603, 2006. 

[57] J. H. Swart. Hopf bifurcation and the stability of nonlinear age-dependent population models. 
Comput. Math. Appl., 15(6-8):555-564, 1988. Hyperbolic partial differential equations. V. 

[58] S. K. Tumuluri. Age- Structured Nonlinear Renewal Equations. PhD thesis, Paris, 2009. 

[59] C. Walker. Prion proliferation with unbounded polymerization rates. In Proceedings of the Sixth 
Mississippi State-UBA Conference on Differential Equations and Computational Simulations, 
volume 15 of Electron. J. Differ. Equ. Conf., pages 387-397, San Marcos, TX, 2007. Southwest 
Texas State Univ. 

[60] G. F. Webb. Theory of nonlinear age-dependent population dynamics, volume 89 of Monographs 
and Textbooks in Pure and Applied Mathematics. Marcel Dekker Inc., New York, 1985. 



35 



